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Abstract. We study a sandpile model on the set of the lattice points in a large lattice polygon. A 
small perturbation i/j of the maximal stable state /j. = 3 is obtained by adding extra grains at several 
points. It appears, that the result of the relaxation of ip coincides with {i almost everywhere; the set 
where ip° ^ is called the deviation locus. The scaling limit of the deviation locus turns out to be a 
distinguished tropical curve passing through the perturbation points. 

Nous considerons le modele du tas de sable sur I’ensemble des points entiers d’un polygone entier. En 
ajoutant des grains de sable en certains points, on obtient une perturbation mineure de la configuration 
stable maximale = 3. Le resultat ip° de la relaxation est presque partout egal a fi. On appelle lieu 
de deviation Tensemble des points oil ip° ^ fi. La limite au sens de la distance de Hausdorff du lieu de 
deviation est une courbe tropicale speciale, qui passe par les points de perturbation. 


1. Introduction 

1.1. Sandpile model on a finite set. Consider the standard lattice on the plane. For t; € Z^ we 
denote by n{v) the set of all four closest points to x; in Z^. Let F be a finite subset of Z^. A state (or a 
configuration) G Z^ of a sandpile on F is a non-negative integer-valued function on F. For a state (j) and 
x; C F we interpret (j){v) as a number of sand grains at the site v. For each x; C F we define the toppling 
operator : Z^ —)■ Z^ at v given by 

Ty(j) = (j) - A6y + y] Syj, 

w^Tr\n{v) 

where 6y is a function on the lattice defined to be 1 at and 0 otherwise. A toppling T„ is called legal for 
a state (j) if (f>(y) > 4, i.e. if Ty(j) is again a state. In this case we think of T„ as a redistribution of sand 
from the overfilled site v to its neighbors. If some neighbors are missing in F, i.e. n{v) ^ F, then at least 
one grain leaves the system after the toppling. If (j){v) < 4 for all n G F, then (j) is called a stable state. 
The state ^ which is defined to be equal to 3 at every point of F is called the maximal stable state. 

A relaxation for a state ^ is a sequence of states 4> = (l> 0 i (I’ll ■ ■ ■ such that (pi+i is the result of 
applying a legal toppling to (pi and (pm is a stable state. It is well known that for any state p there exists a 
relaxation and the last state pm depends only on p (see [Hll]). We denote pm by p° and call it the result 
of the relaxation of p. Informally, in order to find the result of the relaxation it doesn’t matter which legal 
topplings to apply at each step of a relaxation sequence. 

1.2. Motivation. Let II C be a non-degenerate (i.e. of non-zero area) lattice polygon. Let F be 
the intersection of 1? with II. Consider a set P C F. We add extra grains to the state p, at the points 
P = {Pi)P 2 ) • ■ • ,Pn}- After the relaxation, this gives the state ^ = (3 -I- J2peP ^p)°■ 

In the examples shown in Figures |1|2| and we see that the set of points where p is not maximal 
constitutes some sort of a graph passing through P. As it is seen in Figures and the picture becomes 
more regular when the cardinality of P is small with respect to the size of H. In the next section we state 
certain theorems formalizing this concept. Some of this results and their far going generalizations will 
appear in [7]. This short note can be seen as an introduction to the subject and an announcement of these 
results. 
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Figure 1 . Snapshots during the relaxation for the state </> = 3 on a square after adding an 
extra grain at one point p (the big grey point). Black rounds represent x with > 4, 
black squares (which are arranged along the vertical and horizontal edges on the final 
picture) represent the value of sand equal to 2, white rounds (arranged along diagonals 
on the final picture) are 1, and whites cells are 3. Rare cells with zero grains are marked 
as crosses, one can see them during the relaxation on the vertical and horizontal lines 
through p. The value of the final state at p is 3. 

Instantanes pendant la relaxation de la configuration </) = 3 sur un carre, apres ajout 
d’un grain additionnel au point p (le gros point gris). Les ronds noirs representent x avec 
(j){x) > 4, les carres noirs (qui se trouvent sur les aretes verticales et horizontales dans la 
configuration finale) representent les points oii le nombre de grains de sable egale 2. Les 
cercles (sur les diagonales dans la configuration finale) possedent un grain de sable et les 
cellules blanches en ont 3. Les cellules rares avec zero grains sont marquees par des croix; 
on pent les voir pendant la relaxation sur les lignes verticales et horizontales passant par 
p. La valeur au point p dans la configuration finale est 3. 


1.3. Related activities. What we described above is quite close, at least at the basic level, to the original 
Back-Tang-Wiesenfeld model [T]. The fundamental difference is that their framework had a probabilistic 
flavor. Their subject was a stochastic growth of a sandpile simulated by an iterative process of adding an 
extra grain at a random site followed by a relaxation. 

In the paper it has been proven that for any initial configuration most of the states stop occurring 
after some step of the Back-Tang-Wiesenfeld-process and the other states (the recurrent states) occur 
with equal probabilities. Dhar has also shown that the set of all recurrent states forms an abelian group 
(called the sandpile group) under the action (p(B tp = {(f> + ip)° ■ Study of this group is an important and 
fruitful part of the theory of abelian sandpiles [10] . Unfortunately, it is little known about the structure 
of particular elements of the sandpile group [5]. There is a simple characterization of a recurrent state. A 
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stable state (j) is recurrent if and only if it can be represented as (j) = {3 + 1 />)° for some state tp. Therefore, 
the states that we study in this paper are recurrent. 

Our initial inspiration was the definite presence of a tropical curve in Figure 2 (which coincide with 
our Figure]^ in where the authors considered the case where is a square and P consists of just 
one point. The same authors experimentally found more pictures in Section 4.3 of [4] (Figure 3.1 in [12) 1 
for the case of others fl and P; these pictures contain tropical curves too. Note that they considered a 
bit different process: after dropping a grain of sand to a point p and subsequent relaxation, they remove 
a grain from p. After a number of computer simulations, it appeared that the toppling function for the 
relaxation of 3 + J2p(=p is almost piece-wise linear, it was also observed in [Tl], where it was conjectured 
that tropical geometry could be useful in this problem. These articles, mu, containe a lot of observations 
and heuristic arguments. In order to formulate these observations rigorously, we use scaling limits. It 
makes our results particularly close to isiaiisi where the scaling limit of the states {NSo)° was shown to 
exist. Then, the proofs appeared to be rather cumbersome, but evolving into a deep connections between 
sublinear order integer-valued discrete harmonic functions, dynamics of polytopes, and tropical geometry, 
see [7] for details. 


2. The results 

Let n be a non-degenerate lattice polygon and P be a finite non-empty subset of 12°. For any N G N 
consider a set Fat = n Z^. Denote by [p] G 1? the coordinate-wise rounding down of a point p S K^. 
Define the state = (3 -I- '^[Arp])° on Fat and the deviation set 

Cat = -^{v G rAr|<(>iv(ii) < 3}. 

Experimental evidence suggests that when N grows, the shape of sets Cat C D stabilizes, see Figure 

Theorem 1. The sequence of sets Cn C ft Hausdorff converges to a set Cn,p- The set C'n.p is a planar 
graph passing through the points P. Each edge of Cq p is a straight segment with a rational slope. 

Denote by Cn.p the closure of Cq^p n D°. We have a canonical mapping associating the graph Cq^p to 
the configuration of points P. Theorem ^ later represents Cq^p as a solution of some variation problem. 

Eirst of all we note that Cn^p is a weighted graph, i.e. there is a canonical choice of weights for its 
edges. This choice comes from averaging an amount of sand in (pN along edges of Cq^p, Namely, we define 
a sequence of functions —)• K given by (j)'j^{x) = N{3 — (pNiWx])) if [A^a;] G Fat and p'j^ix) = 0 

otherwise. 

Theorem 2. There exists a *-weak limit p of the sequence Moreover, there exists a unique assignment 
of weights rUg T g pj fgj- the edges of Cq,p such that for all smooth functions ip supported on D° 

P{iP) = lim [ (P'j^iP = ^ \\le\\m^'^ [ p), 

where E is the set of all edges of Cn,p and 4 is a primitive vector of e G P, i.e. the coordinates of 4 are 
coprime integers and 4 is parallel to e. 

In other words, edges with high weights are “wider” than edges with low weights (see Figures]^ and [^. 
Theorem motivates the following definition. 

Definition 1 ([T^). Let C be a weighted planar graph whose edges have rational slopes. We define the 
tropical symplectic area of C by 

Area{C) = ^ ||4||™g||e||. 

eeE 

It appears that Ca,p minimizes its area in the class of weighted graphs that we call id-tropical curves. 

Definition 2. A finite planar weighted graph C C il is called an D-tropical curve if 

• each edge e of C has a rational slope and an integer weight mg, 

• the intersection of C with dil is the set of all vertices of D, 
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Figure 2. Triangular domain and two perturbation points. The scale of last picture is 
four times bigger than the scale of the first and two times bigger than the scale of the 
second. 

Un domaine triangulaire, deux point de perturbation. L’echelle de la figure de droite est 
quatre fois celle de la figure de gauche et deux fois celle de la figure centrale. 


• if u S 11° is a vertex of C, then 

( 1 ) ^ rriele = 0 , 

eGB„ 

where Ey is the set of all edges of Cq^p incident to v and le is the primitive vector oi e G Ey 
oriented out of v, 

• there exist a unique labeling dg G Z>o for each side s of fl such that for each vertex u of fl 

(2) ^ ^ — dg^lsi T ds2^S2^ 

eeE^ 

where si and S 2 are the two sides of H incident to v, Ig. is the primitive vector of Si oriented out 
of V. 

Condition 0 is well known in tropical geometry under the name balancing condition, see m- We call 
([^ the outer balancing condition. 

The condition Q implies that an Sl-tropical curve can be always represented as the intersection of fl 
with a tropical curve (see [inE]). A plane tropical curve is a planar graph, whose edges are intervals 
with rational slopes and prescribed positive integer weights. The edges are allowed to be unbounded, 
the number of edges and vertices is finite and at each vertex the balancing condition is satished. If C 
is a subgraph of a tropical curve, then its tropical symplectic area can be seen as the limit of the usual 
symplectic areas for a family of holomorphic curves that degenerates to C (see 13 US])- On the other hand, 
this area can be seen as a particular normalization for the Euclidean length of C. In this interpretation 
C'p,a gives a solution to the analog of the Steiner tree problem. 

Theorem 3. The graph Cq p is an 0-tropical curve. Furtherefore, Cq^p has the minimal tropical sym¬ 
plectic area among all fl-tropical curves passing through the configuration of points P. 

In certain cases, the curve Cn,p can be found by means of this minimization property. This happens 
essentially when 0 is “simple” enough and the number of points |P| is small. In general, for any lattice 
polygon n there are generic configurations P for which Cn^p is not a unique minimizer of the symplec¬ 
tic area. Fortunately, the curve Cq^p can be characterized by the property of minimizing its tropical 
polynomial. 

A tropical polynomial is a function F: > K given by 

F{xi,X 2 )= min (kxi + 1x2 + aki), 

(fc.z)eA 
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Figure 3. The result of adding a grain in two different scales. The polygon 17 is a square 
with edges having slopes 2 and 1/2. 

Resultat de I’ajout d’un grain de sable a deux echelles differentes. Le domaine 17 est un 
carre avec cotes de pente 2 et 1/2 respectivement. 


where A is a finite subset of 1? and ak,i € K. For a tropical polynomial F we consider its corner locus C, 
i.e. the set of points where F is not locally linear. The set (7 has a natural structure of a tropical curve 
m and we say that the curve C is defined by the polynomial F. 

Consider an 17-tropical curve C. Consider an extension of C to a tropical curve C", i.e. C = C" C 17. 
The outer balancing condition (§ implies that there exists a unique tropical polynomial F' vanishing at 
517 which defines C. It is easy to see that F = F'|o does not depend on C. Therefore, we call F the 
ri-tropical polynomial of the curve C. Denote by Fn^p the 17-tropical polynomial of Cq,,p- 

Theorem 4. The polynomial Fq p is the point-wise minimum of all polynomials F of all 17-tropical curves 
that pass through P. 

Theorem can be seen as a manifestation of the least action principle for sandpiles [6] . Consider a 
state (/ on a finite set T C Consider a set % of functions F[: 1? ^ Z>o such that F[ vanishes outside 
r and (j) -I- AH < 3 on T, where AH is the discrete Laplacian of H given by 

AH{v) = -4F(u) -h E «(“-)■ 

w^n{v) 

Denote by F the point-wise minimum of all functions H G H. Then (/)° = </ -|- AF. Furthermore, F(u) 
counts the number of topplings at u € in any relaxation sequence for (j). Therefore, we call the function 
F the toppling function (or odometer) of the state </. 

Let F/v be the toppling function for the relaxation of 3 -I- J2pep on Tat = TV17 C Z^. Consider a 
function F^v: 17 —?► M given by 

Fn{x) = —Fn{[Nx]). 

Theorem 5. The tropical polynomial Fq p is the scaling limit of toppling functions, i.e. the sequence F„ 
converges pointwise to Fq^p. 

This theorem essentially implies all the previous ones. The main idea of the proof of Theorem is 
that the functions F^r are bounded by the piecewise linear function Fn_p and are harmonic on the regions 
where (/jn = h- Then, F^r is harmonic almost everywhere for large N, and its growth is tied by linear 
functions. This finally implies that Fj^ is a piecewise linear function almost everywhere. Detailed proofs 
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and generalizations of these results will appear in the paper [7]. In this forthcoming paper, the lattice 
polygon will be replaced by an arbitrary convex domain. In this case, all the convergence results still 
hold without significant changes. The set Ca,p is still a union of straight edges, but the number of edges 
is usually infinite. It appears that is a tropical analytic curve and is a tropical senes. 
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